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The CoGeNT collaboration has recently published results from a fifteen month data set which 
indicate an annual modulation in the event rate similar to what is expected from weakly interacting 
massive particle interactions. It has been suggested that the CoGeNT modulation may actually 
be caused by other annually modulating phenomena, specifically the flux of atmospheric muons 
underground or the radon level in the laboratory. We have compared the phase of the CoGeNT data 
modulation to that of the concurrent atmospheric muon and radon data collected by the MINOS 
experiment which occupies an adjacent experimental hall in the Soudan Underground Laboratory. 
The results presented are obtained by performing a shape-free \ 2 data-to-data comparison and 
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from a simultaneous fit of the MINOS and CoGeNT data to phase-shifted sinusoidal functions. 
Both tests indicate that the phase of the CoGeNT modulation is inconsistent with the phases of 
the MINOS muon and radon modulations at the 3.0 a level. 



PACS numbers: 95.85.Ry, 29.40.Mc, 95.35.+d 



I. INTRODUCTION 

Numerous astrophysical observations strongly support 
the existence in our galaxy of a cold dark matter 
halo, that may consist of Weakly Interacting Massive 
Particles (WIMPs) 0, [|[ . The principal search mode of 
direct WIMP detection is the identification of an O(keV) 
nuclear recoil produced by WIMP-nucleus elastic scat- 
tering. Since the speed of the Earth relative to the dark 
matter halo varies depending on the Earth's velocity 
with respect to the Sun, the dark matter detection rate 
is expected to demonstrate a annual modulation. This 
modulation is expected to be at a maximum (minimum) 
on June 2 (Dec. 2) with an amplitude between a few 
and 20%, assuming the standard halo model The 
CoGeNT @, 0, DAM A/LIBRA @ and CRESST-II @ 
collaborations have all reported an excess of events 
above all known backgrounds. The CoGeNT and 
DAM A/LIBRA [1, [1(3] collaborations have also claimed 
evidence for annual modulations in their event rates 
at 2.8cr and 8.9 a respectively. Fits to the available 
data favor a light WIMP with mass ~10GeV/c 2 and 
spin-independent cross- section l(T 41 -l(r 39 cm - EMI. 

The null observations by CDMS-II (HUTU, 
XENON100 [IMi and EDELWEISS [20| exclude 
much of the allowed WIMP signal regions mentioned 
above [2l[. The tension between these exclusion 
limits and the positive observations can be signifi- 
cantly reduced, but not removed, when taking into 
account experimental [22. 12311 and astrophysical un- 
certainties [H [TU, [TH, This tension has led 
to suggestions that the CoGeNT and DAM A/LIBRA 
modulations are due to conventional annual phenom- 
ena [111 [29[ . The atmospheric muon rate and the radon 
level in the underground experimental hall modulate 
annually. Signals that can simulate dark matter in- 
teractions may be produced by (a,n) reactions from 
radon decay in the active volume or by nuclear recoils 
from spallation neutrons originating from atmospheric 
muon interactions. The CoGeNT collaboration has 
stated that contamination from these backgrounds is 
small compared to the observed signal 0, [3(|. The 
MINOS experiment monitors both of these quantities in 
an adjacent experimental hall to that of the CoGeNT 
experiment in the Soudan Underground Laboratory. In 
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this paper we compare the modulations of the CoGeNT 
event rate data to that of the atmospheric muon rate 
and radon level data collected at the same time by the 
MINOS experiment. 

The annual modulation of the muon flux deep 
underground has been observed by many different 
experiments (31rl3fjj . The similarities of the amplitudes 
and phases of the modulations observed in the LVD 
muon Q and DAMA/LIBRA data sets motivated 
the hypothesis that modulation in the latter may be 
muon-induced. It has been suggested that spallation 
neutrons or long-lived activated isotopes produced by 
these muons ma y b e responsible for the DAMA/LIBRA 
modulation [28l . |29| | . This now seems unlikely as recent 
detailed comparisons of the DAMA/LIBRA modulation 
to that of the muon fluxes measured by LVD [34[ , Borex- 
ino [33[ and MACRO [3l| , all in the Gran Sasso National 
Laboratory (LNGS), have shown that the phases of 
the two modulations differ significantly |37H39l |. This 
conclusion does not preclude the possibility that the 
CoGeNT modulation, or a significant fraction thereof, is 
due to muon related processes. 

The phase of the modulation of the muon flux can 
vary substantially depending on geographic location 
and calendar year since the flux is strongly correlated 
with the effective atmospheric temperature (3l| - [36| . 
Therefore, to be able to reject with high confidence the 
muon hypothesis as the source of the CoGeNT modu- 
lation, the muon data must be collected concurrently 
with the CoGeNT data and in close proximity to the 
CoGeNT detector. The muon data collected by the 
MINOS experiment fulfill these criteria. Similarly to the 
DAMA/LIBRA muon studies [I?}, we compare the phase 
of the observed MINOS muon modulation to that of the 
CoGeNT data modulation. Comparisons of the CoGeNT 
data to non-concurrent MINOS muon data [36| . and 
indirectly to effective temperature variations, have been 
presented in Ref. [39| and indicate that the data sets are 
not correlated. 

We note that the 16.6% amplitude of the CoGeNT 
event rate modulation Q is significantly larger than 
the ^2% amplitude of the MINOS muon rate mod- 
ulation (36[. This difference suggests that the muon 
temporal variation cannot fully account for the observed 
CoGeNT modulation. In this paper we examine the 
relative phases of the two modulations which provides 
an independent test of the potential correlation between 
the CoGeNT and MINOS muon data sets. 
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The radon level in the Soudan Underground Labora- 
tory is at a maximum (minimum) in the summer (winter) 
months due to the pressure gradients created by the 
relative temperature differences between the air in the 
laboratory and that on the surface [4(|. In the MINOS 
cavern we have observed that the radon concentration 
varies by a factor of six over the year, corresponding to 
a modulation amplitude of ^60%. A large modulation 
amplitude could therefore be introduced into the Co- 
GcNT data by even a small amount of contamination 
from this background. 

The radon progeny also modulate with a one year 
period T, but do so with a delayed phase and reduced 
amplitude. The decays between 222 Rn and 210 Pb occur 
very quickly (^minutes) and therefore have negligible 
impact on either the phase or the amplitude. Since 
210 Pb has a half-life of 22 years, its decay and the decays 
of its progenies will not contribute to the modulation. 



has published its results using data collected over a pe- 
riod of 458 days between Dec. 4, 2009 and Mar. 6, 2011 
with a total of 442 live days @ . The data were presented 
in fifteen 30-day intervals and one 8-day interval, then fit 
to a modulation hypothesis of the form: 



R = R 1 + A ■ cos 



2tt , 

jr(*-*o) 



(1) 



where i?o is the mean rate, A is the modulation am- 
plitude and T is the period. The time t is the number 
of days since Jan. 1, 2010. The phase t is the day 
at which the signal is at a maximum. The published 
CoGeNT best fit results are given in the last line of 
Table |IJ The modulation hypothesis is preferred over the 
null hypothesis at 2.8 a. The CoGeNT collaboration has 
released the background-subtracted data set used in this 
analysis to the public. The results of our \ 2 fit of the 
CoGeNT data to Eq. {TJ, discussed further in Sec. JIII A| 
are in good agreement with the published results 



The following Section of the paper discusses the se- 
lected experimental data sets. In Section IlIII we present 
the best fit modulation parameters determined for each 
of these data sets. We then describe the measurements of 
the phase differences between the CoGeNT and MINOS 
muon and radon data sets obtained from a simultaneous 
fit of the data to phase-shifted sinusoidal functions, a 
shape-free x 2 data-to-data comparison and a bin by bin 
correlation test. Section [IVI summarizes our conclusions. 



II. THE SELECTED DATA 

The CoGeNT dark matter experiment [30|, EJ and 
the Far Detector of the MINOS long baseline neutrino 
experiment [42| are located 705 m underground in two 
different caverns of the Soudan Underground Laboratory. 
The MINOS cavern, which houses the MINOS detector, 
is 82 m long, 15 m wide and 13 m high and is oriented 
along the direction of the NuMI neutrino beam (43| . 
The CoGeNT and CDMS-II dark matter experiments 
are located in the Soudan 2 cavern which is similar in 
shape to the MINOS cavern but is 70 m long and is 
oriented north-south. The two experimental caverns are 
connected by an east-west passage on their north side 
and are served by a common ventilation system which 
replaces the lab air several times per hour. 



A. The CoGeNT Data 

CoGeNT is an experiment for direct detection of dark 
matter which employs a 0.44 kg p-type point contact ger- 
manium detector 00, 0- The CoGeNT collaboration 



B. The MINOS Data 

The MINOS Far Detector has been collec ting atmo- 
spheric muon data since August 2003 [12, l45| . The 
experiment also records the radon level in the labora- 
tory air. The MINOS muon and radon data used in 
this analysis were collected between June 4, 2009 and 
Sept. 6, 2011. This collection window is 12 months 
longer than the CoGeNT run period, from Dec. 4, 2009 
to March 6, 2011, allowing the data-to-data comparisons 
described in Sections IIII Bl and MI CI 

The event selection and data quality requirements 
used in this analysis are identical to those in the pre- 
vious study of seasonal muon intensity variation at the 
MINOS Far Detector, with the additional requirement 
that the reconstructed muon track be downward going. 
Restricting the data set to contain only days with greater 
than 10,000 s of live time yields a total of 738 good days 
of atmospheric muon data. These good days include 
449 days which occurred between Dec. 4, 2009 and 
March 6, 2011 inclusive. 

The radon level in the MINOS cavern air, inferred from 
counting the number of alpha decays, is measured every 
hour by a Model 1027 Sun Nuclear Corporation radon 
monitor [46j . A daily measure of the radon level is deter- 
mined by averaging the 24 measurements taken through- 
out the day. The standard deviation of these measure- 
ments, a, is taken to be the error on the daily radon mea- 
surement. While larger than the standard error on the 
mean value, a/y/^A, this choice is more consistent with 
the published accuracy of the radon monitor [46| . There 
are 786 good days during which the radon monitor oper- 
ated continuously throughout the day. These good days 
include 458 days which occurred between Dec. 4, 2009 
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and March 6, 2011 inclusive. The radon monitor was 
moved to different locations in the Soudan Underground 
Laboratory and cross calibrated with other detectors run- 
ning simultaneously. This demonstrated that the radon 
level does not vary spatially in the laboratory to within 
the resolution of the monitor. Thus the radon levels 
measured in the MINOS cavern can be used to evalu- 
ate whether the CoGeNT data are correlated with the 
radon level in the Soudan cavern. 
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cantly different than that of the MINOS muon or radon 
data we can infer that they are likely not causally related. 

The most common approach in the literature to 
evaluating potential correlations, and discussed here in 
Sec. IIII A[ is to fit the data to Eq. (Q} and compare 
the phases and periods of the best fits. The CoGeNT 
and DAMA/LIBRA modulations are a good fit to a 
cosine function. This is the expected signature for an 
isothermal dark matter halo. The true form of the 
modulation may be more complex as it is dependent on 
assumptions made regarding the velo city distribution of 
the dark matter particles in the halo [HI, HtJ . The muon 
modulation is not fit well by a cosine function [HI, H^] . 
The muon and radon modulations are correlated with 
atmospheric temperatures. Therefore, their modulations 
are cyclical but not necessarily sinusoidal. Imposing 
such constraints onto the data may bias the results of 
the cosine based fit comparison. We address this concern 
in Sections IIII Bl and IIII CI by performing shape-free 
data-to-data comparisons that allow us to evaluate the 
phase differences and potential correlations regardless of 
the underlying functional forms of the modulations. 



A. Cosine \ 2 Test 



FIG. 1: The residuals of the MINOS Far Detector muon rates, 
radon levels and CoGeNT event rates as a function of time. 
The MINOS muon and radon data have been scaled by factors 
of 10 and one-half respectively to fit on the same graph and, 
for this figure, use the same binning as the CoGeNT results. 
The vertical dashed lines indicate the start of a new calendar 
year. The arrow marks the date where a dark matter signal 
is expected to peak. 

The MINOS muon rate and radon level residuals, 
and the CoGeNT event rate residuals, arc plotted 
as a function of time in Fig. [TJ The CoGeNT event 
rate residuals are calculated with respect to a mean 
rate of 97.7 events/30 days. The MINOS muon rate 
residuals are calculated with respect to a mean rate of 
(0.4431 ± 0.0001) Hz. The MINOS radon level residuals 
are calculated with respect to a mean level (11.94 ± 
0.11)pCi/l. All three data sets possess clear modulation 
signatures. In the following section we quantify any 
potential correlations between these modulations. 



III. MODULATION COMPARISONS 

If the CoGeNT modulation is caused by either the 
muon or radon backgrounds then it should modulate 
with the same shape as those backgrounds. Therefore, 
if the phase of the CoGeNT modulation is signifi- 



The nominal modulation parameters for the CoGeNT 
and MINOS muon and radon data sets were determined 
by performing a x 2 fit test of Eq. JT]) to the data 
described in Sec. [TT| and shown in Fig. [T] The results of 
these fits are given in Table HI The confidence limit con- 
tours for the best fit phase and period are shown in Fig. [2] 

Our fit to the CoGeNT data is in good agreement 
with the published results [f| and disfavors the null 
modulation hypothesis at 3.1 a. The significance with 
which we exclude the null modulation hypothesis is 
defined as the square root of the difference between 
the x 2 value of the best fit point and that of the null 
modulation hypothesis. This definition is different from 
that used in the published CoGeNT analysis and gives 
a slightly stronger exclusion. The small differences 
between our best fit values to the CoGeNT data and the 
published CoGeNT best fit values may be explained by 
the assumption in our fits that the CoGeNT errors are 
uncorrclated. 

The two apparent occurrences of sudden stratospheric 
warming events [48| in early 2010 and early 2011, which 
temporarily increased the muon rate, drive the large 
X 2 for the muon fit and cause the best fit period to 
be significantly smaller than one year. If the complete 
MINOS muon data set, August 2003 to April 2012, 
is fit, minimizing the impact of short term fluctua- 
tions, a period much closer to one year is obtained, 
T= (364.5 ± 0.3) days, and the phase remains unchanged. 
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Data 




X 2 / N -d-o.f. Mean Rate 


Amplitude Period 


Phase 


Date of 






[•Ho] 


[A,%] [T.days] 


[to, days] 


Maximum 


Best fit modulation parameters assuming a fixed period of 365.25 days. 


Muon 




1909 / (449-3) (0.4428 ± 0.0001) Hz 


1.25 ± 0.03 365.25 


182.8 ± 1.7 


July 1 


Radon 




176 / (458-3) (11.9 ± 0.1) pCi/1 


57.7 ± 0.9 365.25 


215.0 ± 1.1 


Aug. 3 


CoGeNT (Our 


Fit) 


6.6 / (16-3) (97.9 ± 3.6) counts/30 days 


16.9 ± 5.4 365.25 


108.4 ± 16.9 


Apr. 18 






Best fit modulation parameters without 


a fixed period assumption. 




Muon 




1788 / (449-4) (0.4431 ± 0.0001) Hz 


1.37 ± 0.04 317.2 ± 3.2 


187.3 ± 1.4 


July 6 


Radon 




176 / (458-4) (12.0 ± 0.1) pCi/1 


57.7 ± 0.9 367.4 ± 3.5 


215.2 ± 1.1 


Aug. 3 


CoGeNT (Our 


Fit) 


6.4 / (16-4) (97.7 ± 3.6) counts/30 days 


16.7 ± 5.4 348 ± 42 


113.7 ± 17.9 


Apr. 23 


Published CoGeNT modulation parameters [6]. 


CoGeNT 




7.8 / (16-4) N/A 


16.6 ± 3.8 347 ± 29 


~115 ± 12 


Apr. 25 



TABLE I: The best fit results produced by fitting the MINOS muon rate, radon level and CoGeNT event rate data to Eq. {TJ. 
The fits reported in the first three rows of the table have been performed with the period fixed to lyear (365.25 days). The 
last column gives the dates in 2010 at which the fits to the data are at a maximum. 
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FIG. 2: Confidence limit contours for the period and phase as 
determined by fitting the CoGeNT event rate and MINOS Far 
Detector muon rate and radon level data to Eq. {TJ. The best 
fit values are given in Table [i] The horizontal and vertical 
black lines mark the expected period and phase for a dark 
matter signal. 



The best fit phase differences Sto between the CoGeNT 
phase and the MINOS muon and radon phases are deter- 
mined by minimizing: 



number of live days concurrent with the CoGeNT data 
collection period. The second term is the contribution 
from the CoGeNT event rate data. R b,M,i {Rob,c,i) is 
the i th observed MINOS (CoGeNT) data point.' a M ,i 
and ere. j are the uncertainties on the MINOS and 
CoGeNT data points respectively. R ex is the expected 
value, as determined by Eq. ([T]), assuming the given 
modulation parameters and Sto is defined as the phase 
of the CoGeNT data minus the phase of the MINOS 
data. The x 2 , as a function of this phase difference, is 
determined by minimizing the x 2 over the MINOS mean 
value Ro,m, the amplitude Am and phase to and the 
CoGeNT mean value i?o,c, amplitude Ac', and, for some 
fits, a common period T. 

Figure [3] shows the A% 2 curves, as a function of Sto, for 
the simultaneous fits of the MINOS and CoGeNT data 
to Eq. ([2]) assuming a common period of one year. The 
best fit phase differences are (-75 ± 18) days and (-110 ± 
18) days for the comparison to the muon and radon data 
respectively and (-67 ± 17) days and (-112 ± 18) days 
respectively when minimizing the x 2 over the period T. 
The statistical significance at which equivalent phases for 
the MINOS and CoGeNT data can be excluded is given 
by the square root of the Ax 2 difference between the best 
fit point and the value at Sto = 0. As can be seen from 
Fig. [3] the phases of the MINOS muon and radon data 
are inconsistent with the phase of the CoGeNT data at 
3.0 a and 3.1 a respectively. 



2/r, \ (Rob.M.i — Rex(Ro,M, Am, to, T)) 2 

X (Sto) = > ; 2" J ( 2 ) B. Shape-Free x Test 

16 (Rob,c,i - Rex(Ro,c A c,to + St a ,T)) 2 In this section we determine the relative phase St be- 

0-2, _ tween the MINOS and CoGeNT data sets, without an a 

1-1 priori assumption regarding their shape, by calculating 

The first term in Eq. (J5|) is the x 2 contribution from the the x 2 difference between their respective modulations. 

MINOS muon rate or radon level data where Nm is the The x 2 difference, assuming a common binning, is defined 
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FIG. 3: The Ax 2 distributions comparing the phases of the 
MINOS muon rate and radon level data to the phase of the 
CoGeNT event rate data using Eq. ([2]). The A^ 2 curves are 
calculated with respect to their x 2 minima. The flattening of 
the Ax 2 curves indicate that these exclusions are limited by 
the confidence with which the CoGeNT data can exclude the 
null modulation hypothesis. 



FIG. 4: The Ax 2 distributions comparing the phases of the 
MINOS muon rate and radon level data to the phase of the 
CoGeNT event rate data using Eq. @. The Ax 2 curves are 
calculated with respect to their x 2 minima. 



C. Correlation Test 



as: 

2rx+ i (R c ,i - f ■ RmA^o)) / q \ 
X (Sto) = 2^ 2 2 ■ (3) 

Rm,i (Rc,i) is the i th MINOS (CoGeNT) residual and 
(?m i and ac i are the uncertainties on the MINOS and 
CoGeNT residuals respectively. We marginalize over the 
difference in amplitudes, for each Sto, by minimizing the 
X 2 over a positive definite multiplicative factor /. If the 
data have similar underlying forms, we expect the \ 2 to 
be a minimum when the phase difference between them 
is zero. The x 2 values, as a function of St, are deter- 
mined by shifting the time-axis of the MINOS data by 5t 
days and recalculating Eq. Figure 0] shows the Ax 2 
curves as a function of the MINOS data offset, which is 
equivalent to the relative phase 5to. The curves are not 
smooth due to statistical fluctuations in the data. By off- 
setting the MINOS data we vary the number of MINOS 
live days which overlap the CoGeNT data. To ensure 
that each subset of MINOS data, for every St, contains 
the same number of live days we substitute the historical 
daily average of that date for those days which do not 
pass the live-time selection criteria. The best fit phase 
differences between the CoGeNT data and the MINOS 
muon (radon) data, corresponding to the minimum of the 
Ax 2 curves in Fig. H are -83± 25 days (-123±j|days). 
The statistical significance, as defined in Sec. IIII A[ at 
which equivalent phases for the CoGeNT and MINOS 
muon (radon) data are excluded is 2.9c (3.0 a). 



Residual muon or radon backgrounds in the CoGeNT 
data could cause a correlation between the CoGeNT 
modulation and the MINOS muon and radon modula- 
tion measurements. The degree of correlation has been 
evaluated using Pearson's coefficient of correlation, cal- 
culated as: 



1 (Rob,M,i — Rob,Al)(Rob,C,i — Rob,c) 

Nc — 1 CFMOC 



where Nc is the number of bins, R b,M and R b,c are the 
average values of the MINOS and CoGeNT data sets, gm 
and ac are the standard deviations of the points compris- 
ing the MINOS and CoGeNT data sets respectively. The 
correlation coefficients, and their Fisher transforms (49| . 
are given in Table HU 



Data Set 


Correlation 


Fisher 




Coefficient (p) 


Transform 


CoGeNT vs Muon Data 


0.19 


0.19 ± 0.28 


CoGeNT vs Radon Data 


-0.29 


-0.30 ±0.28 



TABLE II: The coefficients of correlation, and their Fisher 
transforms, calculated between the CoGeNT event rate data 
and the MINOS muon rate and radon level data. Both data 
sets being compared are consistent with no correlation at ~1<t. 

Even if there is no causal relationship between the 
observed MINOS muon and radon modulations and 
the CoGeNT modulation, there will be some corre- 
lation between these data sets as they all follow an 
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approximate sinusoidal variation. The expected value of 
the correlation is related to their relative phases. For 
example, if the phase difference between two periodic 
data sets is smaller (larger) than one-quarter of the 
period, the correlation should be positive (negative). 
One can therefore infer from the results in Table [TT] 
that the effective phase difference between CoGcNT 
and the MINOS muon data is near to but less than 
365.25/4 days, while between CoGeNT and the MINOS 
radon data it is near to but more than 365.25/4days. 

To verify whether the calculated correlations are 
consistent with the observed modulation phases we 
generated a series of pseudo-experiments. Sampling 
from two cosine curves, with the precision and binning 
of the CoGeNT and MINOS data sets and amplitudes 
taken from Table HI we calculated the Fisher transform 
as a function of the phase difference between the two 
curves. We find that the observed values of the Fisher 
transforms in Table [TO (0.19 ± 0.28) and (-0.30 ± 
0.28), correspond to phase differences of -77l47days and 
-117137 days respectively. These values are consistent 
with the phase differences calculated in the preceding 
sections. 



IV. CONCLUSION 

We have performed a comparison of the modulation 
phases observed in the CoGeNT and MINOS atmo- 
spheric muon and radon data, all collected concurrently 
between Dec. 4, 2009 and March 6, 2011 in the Soudan 



Underground Laboratory. We have presented the 
results of a shape-free data-to-data comparison which 
indicate that the phases of the CoGeNT data and the 
atmospheric muon and radon data are different by 
-83tfdays (2.9 a) and -123±^days (3.0 a) respec- 
tively. The calculated correlation coefficients between 
the CoGeNT and MINOS data sets are statistically 
consistent with the no-correlation hypothesis. The 
cosine fit test measures the phase difference between 
the CoGeNT and MINOS muon data sets to be (-75 ± 
18) days, inconsistent at 3.0 a, and between the CoGeNT 
and MINOS radon data sets to be (-110 ± 18) days, 
inconsistent at 3.1a. The similarity between the results 
of both these tests indicate that no significant bias is 
introduced when imposing a sinusoidal shape on the 
data. It is also clear that our exclusions are limited 
by the degree to which the CoGeNT data exclude 
the null modulation hypothesis. Based on the studies 
described above, it appears unlikely that muon or radon 
related processes contribute significantly to the observed 
CoGeNT modulation. 
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